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We show that in a general = 1 supergravity with iV 3> 1 scalar fields, an exponentially small 
fraction of the de Sitter critical points are metastable vacua. Taking the superpotential and Kahler 
potential to be random functions, we construct a random matrix model for the Hessian matrix, which 
is well-approximated by the sum of a Wigner matrix and two Wishart matrices. We compute the 
eigenvalue spectrum analytically from the free convolution of the constituent spectra and find that in 
typical configurations, a significant fraction of the eigenvalues are negative. Building on the Tracy- 
Widom law governing fluctuations of extreme eigenvalues, we determine the probability P of a large 
fluctuation in which all the eigenvalues become positive. Strong eigenvalue repulsion makes this 
extremely unlikely: we find P oc exp(— cA'^^), with c, p being constants. For generic critical points 
we find p ^ 1.5, while for approximately-supersymmetric critical points, p ^ 1.3. Our results have 
significant implications for the counting of de Sitter vacua in string theory, but the number of vacua 
remains vast. 
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1 Introduction 



Perhaps the most pressing question in string theory is whether the theory admits solutions 
consistent with all observations. In light of the discovery of the acceleration of the universe, 
it is essential to pursue de Sitter solutions of string theory, and to understand whether these 
solutions are so numerous that they can account for the smallness of the vacuum energy. After 
a decade marked by significant advances in understanding flux compactifications [T], there is 
now compelling, but still largely indirect, evidence for the existence of a vast landscape of 
metastable de Sitter vacua. Direct enumeration of explicit de Sitter vacua remains a distant 
goal. 

The cardinal difficulty in constructing de Sitter solutions is that in the absence of super- 
symmetry, the scalar potential can have instabilities along one or more directions in the scalar 
field space. When the number of fields is large - which is both generic in Calabi-Yau fiux 
compactifications, and indispensable for providing an astronomical number of vacua - direct 
examination of the Hessian matrix of the scalar potential becomes impractical. This impasse 
motivates a statistical approach, in which the compactification data are taken to be random 
variables. 

As a metastable vacuum is a critical point of the scalar potential at which the Hessian 
matrix is positive definite, it is natural first to ascertain the statistical properties of general 
critical points, and then to characterize the subset of critical points that are in fact metastable 
vacua. In the seminal work [2] , Denef and Douglas formulated this problem for a general four- 
dimensional M = 1 supergravity theory, taking the superpotential W and Kahler potential K 
to be random functions, in a precise sense that we shall review. Denef and Douglas studied the 
eigenvalues of the Hessian matrix % and argued that a significant fraction of critical points are 
metastable vacua. 

In this work we reexamine the stability of de Sitter critical points in a general four- 
dimensional M = 1 supergravity. Our tool is random matrix theory: is a large matrix, 
and a great deal can be said about its eigenvalue spectrum. Moreover, given an ensemble of 
random matrices that typically have negative eigenvalues, the Tracy- Widom theory of fiuctua- 
tions of extreme eigenvalues allows one to compute the probability of drawing a positive-definite 
matrix from the ensemble The key phenomenon is eigenvalue repulsion: a large fiuctua- 
tion through which all eigenvalues become positive generally requires an increase in the local 
eigenvalue density, which is statistically costly. 

We obtain results that depend on the relative sizes of the soft supersymmetry-breaking 
masses and the supersymmetric masses. At a generic critical point, the supersymmetric masses 
are not hierarchically larger than the soft masses, and supersymmetry provides negligible pro- 
tection from instability. We develop a random matrix model for "H and obtain an analytic 
expression for its eigenvalue spectrum. The spectrum has considerable support at negative 
values, so that tachyons are generic. Building on extensions of the Tracy- Widom theory due to 
Dean and Majumdar [H |3], we then argue that the probability P of a fiuctuation rendering % 
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positive-definite is 

P (xexp{-cNP) , (1.1) 

at leading order in large A^, where N is the number of complex scalar fields, and c, p are 
constants, with p ^ 1.5. This is qualitatively similar to the original results of Aazami and 
Easther [6], who obtained p ^ 2 in a simpler model of the Hessian matrix. We conclude that 
an exceedingly small fraction of generic critical points are metastable. 

The more promising regime, as stressed by Denef and Douglas, is that in which approxi- 
mate supersymmetry protects against instabilities. When the soft masses are small compared 
to the supersymmetric masses, the only significant risk of an instability arises from the di- 
rection parameterized by the scalar superpartner of the Goldstinoj^ and we show that at a 
critical point that is generic apart from this requirement of approximate supersymmetry, there 
are almost always two tachyons. We identify a negative-definite contribution to these eigen- 
values that is a manifestation of eigenvalue repulsion between the Goldstino directions and the 
supersymmetrically-stabilized scalars. The eigenvalue repulsion contribution is dominant at 
large N, so that it is extremely improbable that both eigenvalues fluctuate to become positive: 
the probability of positivity again takes the form (1.1), but now with p ^ 1.3. 

Although our results clearly demonstrate that an exponentially small fraction of the critical 
points of a generic random supergravity theory are metastable, this finding in no way precludes 
the existence of a landscape of metastable de Sitter vacua. First of all, the actual number of 
metastable vacua can be extremely large - and in particular, larger than 10^^° - while still being 
exponentially small compared to the number of critical points. Second, our analysis applies 
when K and W are generic functions of all of their arguments. A well-motivated configuration 
violating this assumption is a theory involving two decoupled sectors, containing A'^^^ and 
scalar fields, respectively. If the Nh fields receive large supersymmetric masses and dynamical 
supersymmetry breaking occurs at a lower scale in the sector containing Nl fields, then the 
dominant factor in the number of critical points can be exponential in Nh, while the fraction 
of critical points that are unstable is proportional to exp{—cNf^). Thus, for Nh ^ N^ our 
findings yield only a modest reduction in the number of vacua. 

The organization of this paper is as follows. In ^ we set our notation and review the 
structure of the Hessian matrix "H at a critical point in a general four-dimensional A/" = 1 
supergravity, following [2]. In ^we introduce the ideas from random matrix theory that are 
essential in this work, reviewing the relevant ensembles and assembling results from the theory 
of fluctuations of extreme eigenvalues. In Q we apply these methods to study 7/ at a generic 
critical point. We compute the eigenvalue spectrum analytically and obtain the probability of 
a large fluctuation that renders H positive-definite. In ^ we study T-L at an approximately- 
supersymmetric critical point. We show that eigenvalue repulsion typically leads to two negative 
eigenvalues associated to the Goldstino, and we again compute the probability of a fluctuation 
to positivity. In ^ we discuss extensions of our assumptions, and we illustrate our results in 
the KKLT scenario. We conclude in ^ 



^For extensive investigations of this unstable direction, see [?]■ 
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2 Critical Points in Af = 1 Supergravity 



In this section we will discuss the form of the critical point equation, and describe the structure 

we 



of at a critical point, in a general four-dimensional JV = 1 supergravity. In [2A_ and ^ 2.2 
closely follow the work of Denef and Douglas [2] , reviewing how the critical point equation can 
be written as an eigenvalue equation, and we establish notation for the different contributions 



to the Hessian matrix. Then, in §2.3[ we give a precise definition of a random supergravity, 
whose critical points will be the object of study in ^ and [|5j 

The F-term potential in an A/" = 1 supergravity with chiral superfields is given by 

V = {FaF" -3\W\'^) , (2.1) 

with a = 1, ... A^, in units in which SirG = M'^ = 1. Here = K^'^Fh = K'^^Dj^, where Da 
is the Kahler covariant derivative, DaW = daW + KaW, and K^i is the Kahler metric. We will 
use the shorthand FaF"' = F'^ . 

We consider a set of critical points {g}, satisfying 

daVU = e''(Va{F^)F' - 2FaW'^ = . (2.2) 

Here and henceforth Va denotes the appropriate Kahler and geometrically covariant derivative. 
At any given point q the scalar potential can be simplified by specifying the Kahler gauge such 
that {K)q = and performing an appropriate coordinate transformation such that K^ilq = 5^1- 



2.1 Matrix form of the critical point equation 



The critical point equation (2.2) can be written as an eigenvalue equation of a particular 
Hermitian matrix, A^, formed from the second covariant derivatives of the superpotential [2]. 
Defining 

Zab = T^aFb , (2.3) 



and with = Arg(iy), equation (2.2) can be expressed as 

MF = 2\W\F, 

where 

M = 







ah 



and the 2A^-dimensional vector F is given by 







(2.4) 
(2.5) 

(2.6) 



The eigenvalues of Ai come in real pairs with opposite signs, ±Aa, with Aa > and a = 1, . . . , A^. 
Thus, at any critical point, A4 must have an eigenvalue equal to 2|iy|, and the vector F must 
be proportional to the corresponding eigenvector. In ^3.1.1 we will discuss the spectrum of 
eigenvalues of A/1. 
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2.2 Structure of the Hessian matrix 

For a critical point g to be a metastable vacuum, the eigenvalues of the Hessian matrix of the 
scalar potential evaluated at q must all be positive. Denoting Uabc = T^aPhFc , and writing 
the curvature of the field spac^ in terms of the partial derivatives of the Kahler potential as 
Rabcd = ^af ^c^id ~ ^abcd ~ ^ac^^bde: bosouic mass matrix can be written 

dlV = UabcF'-WZ^b, (2.7) 
dliX = 5,i{F^-2\W\'')-Fah-RalcdF'F^+Z: Z-,,, (2.8) 

where indices are raised and lowered with 5°"^ . The Hessian matrix % is thus given by 

Za"" ^bc - FaFl - R^lc^FT'^ UabcF" - ZabW _ \ ^ 

U abcF'' - Z^iW Z-^ Zbc- FbF-a- R^jj^^^F'^F'^ ) 

+ \[f'^ - 2\wf^ , (2.10) 

where 1 denotes the IN x 2A^ unit matrix. The Hessian matrix is most conveniently analyzed 
in a 'Goldstino' basis in which Fa = d^Fe^"^^. In this basis the critical point equation (2.4) can 
be written 

Zn = 2\W\ 6^(2''^-''^) , Zia' = , (2.11) 

while the components Za>b' remain unconstrained for a',b' = 2, . . . , N. The Hessian matrix can 
be decomposed into constituent matrices as follows: 



where 



"^susy - ( " Q '''' Z " Z ^ ' {2.13) 



be 



- ( ~"r -L ) • <^-^^' 

Hshift = i(f2-2|W^|2) -FV^^I-FV^,!- (2.17) 

A few remarks are appropriate at this point. First, "Hsusy and "H^cs) are positive semidefinite. 
Second, the mass scale mgusy of the supersymmetric masses is set by the eigenvalues of Hsusy, 



^Our sign convention differs from that of [2]. 
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and can be larger or smaller than the scale F that determines the soft supersymmetry-breaking 
masses. The ratio F/nisusy (recall that we have set Mpi = 1) has a significant effect on stability. 
In we will study generic critical points, at which|^ F ~ Tris^sy, and in ^ we will consider 
'approximately-supersymmetric' critical points at which F <^ J^j-susy 

In ^ we will explain how the constituent matrices of "H given in equations (2.13)-(2.17) can 
be identified as — or well-approximated by — elements of classical random matrix ensembles 
with well-known emergent eigenvalue spectra at large N. The distribution of the eigenvalues 
of the Hessian matrix can then be obtained as the free convolution [H] of the eigenvalue distri- 
butions of the constituent matrices, just as the distribution of a scalar random variable that is 
the sum of terms with known distributions can be obtained by the ordinary convolution of the 
constituent probability density functions. 



2.3 Defining a random supergravity 

To proceed further, we must specify the statistical properties of the entries of the matrices 



(2.13)-(2.17) constituting the Hessian matrix Ti. Our fundamental assumption — consistent 
with that of [2j — is that the components of tensors formed by covariant differentiation of 
W and K are independent, identically distributed (i.i.d.) variables drawn from some statistical 
distribution Q. We will occasionally abbreviate this by saying that W and K are random 
functions. Note that this assumption is quite different from taking the entries of "H itself to be 
i.i.d. variables drawn from a distribution Q, which omits the structure and correlations implicit 



in (2.12) 



We will use a) to denote a complex]^ distribution whose magnitude has mean jj and 



standard deviation a, with a uniform distribution for the phase. In ^6.1 we explain that as a 
consequence of the well-known phenomenon of universality in random matrix theory, the precise 
choice of Q is immaterial, provided that the higher moments of Q are appropriately bounded. 



2.3.1 The Kahler potential and its derivatives 

Suppose we took the Kahler potential to be a random function such that in a generic coordinate 
basis, 

KJ,enio,i). (2.18) 

Performing an orthogonal rotation to diagonalize Kg],\g, the resulting eigenvalues will generically 
be of order (see ^for details), and the GL(N, C) transformation required to achieve i^'aj = Sai 
involves rescaling by factors of order A^. To avoid performing this rescaling in all terms involving 
K, we find it convenient to take 

Kj^em^)- (2.19) 

■^See Appendix A for a discussion of the distribution of F/rrisusy in the set of all critical points. 
^The diagonal elements of Hermitian matrices will of course be real. 
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Then, the GL{N, C) transformation leading to K^j, = 5^1 does not involve any A^-dependent 
rescalings. More generally, the choice f2(0, ^) is convenient because the random matrix eigen- 
value spectra presented in ^ then have support in the same domain for all N . 

Next, we need to specify the properties of K^i^ and K^i^g^ at g, in the basis in which 
K^l = Sab- Ideally, the statistics of these objects would follow from a theory of general Kahler 
metrics (see e.g. [9] for a very recent discussion of related issues), but for our purposes it will 
suffice to stipulate that the Kahler potential K is a random function of its arguments, in the 
sense described above. Then, i^^^^ and -ft'afccJ do not take a special form in the basis in which 



the metric is diagonalized, and imposing the critical point equation (2.4) does not change this 
situation. In the 'Goldstino' basis in which Fa = Sj^F e^"^^ at g, we have 



K. 



ab\q 



K. 



ahl\q 



K. 



abll\q 



G 
G 



fi(0, 
1](0, 



(2.20) 
(2.21) 
(2.22) 



The assumptions (2.21), (2. 22) are well-motivated for general Kahler manifolds, but we note 



that there are interesting exceptions, including the special geometry of the vector multiplet 
moduli space m M = 2 supergravity, for which the curvature tensor is given by 



R. 



abed 



K. 



abed 



bde 



^nh^ed + ^ad^eb 



2K TT-pq -p XL _ 
^ -J aep'' bdq i 



(2.23) 



where J-" is the prepotential. Repeating the analyses of §^4][5 with the special geometry rela- 



tionship (2.23) is straightforward, and we find a decreased likelihood of stability compared to 



the more general assumptions (2.21), (2. 22) that are used throughout this work. 



2.3.2 The superpotential and its derivatives 

Turning now to the superpotential and its derivatives, we begin with a warmup in global 
supersymmetry. Fixing a point q in field space and working with canonically-normalized fields 
(pA, A = 1, ... A^, we may write W in the form (momentarily restoring factors of the Planck 
mass for clarity) 

W = M^ w{<f)i/Mpi, . . . , (j)N/Mpi) , (2.24) 

where M is a mass scale. Our assumption is that w is a random function of its dimensionless 
arguments xa = 4>a/Mpi, so that 



dw 



dx, 



G n(/x, a) 



G f2(/i, a) 



(2.25) 



dxAdxB 

etc. At a typical point the various derivatives of the superpotential will be of the same order 
of magnitude, set by the physical effect responsible for the superpotential. (For example, in 
type IIB flux compactifications, the sizes of the superpotential and its derivatives are set by 
the flux scale.) However, atypical points will play an important role in ^ it can happen that 
the superpotential, as well as its flrst derivative F^, are small compared to higher derivatives. 
The result is a signiflcant change in the stability criteria [2]. 
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In supergravity, the relevant expansion around q is in (Kahler and geometrically) covariant 
derivatives of the superpotential. We take 



(2.26) 
(2.27) 
(2.28) 

We will assume that supersymmetry is spontaneously broken by an F-term; for a discussion 



Fa = VaW 


e 


'^susy 


n{o, 


1 


Zab = VaV.W 


e 


^susy 


fi(0, 


1 


Uabc = vjy^v.w 


e 


^susy 


fi(0, 


1 



of the possible effects of D-term energy (cf. |2j), see [6.3 The requirement of nonnegative 
vacuum energy then becomes 

F > V3\W\ , (2.29) 
so that in particular, < 0{F). It is useful to define 



3\W\ 



(2.30) 



so that V = F^(l — u"^), and equation (2.29) translates to w < 1. The stability properties of 
critical points depend on u, so we will repeat our analysis for a collection of fixed values of 
oj e [0, 1]. Thus, \W\ is determined in terms of F and u. For more details on this point, see 
Appendix lAl 



We can now construct a model of a random supergravity by drawing Fa, Zab, Uabc, K 

.(4) 



(3) 
abc'' 



and K^i^^i independently from the distributions specified above, at each point q. As we are 



primarily interested in critical points, we should study the set of points {q} subject to the 
critical point equation (2.4). Such points are not completely generic: equation ( |2.4 ) enforces a 
particular correlation between Zab, W , and Fa, as reviewed in the discussion following equation 



(2.4). Following we will carefully incorporate the restriction implied by equation (2.4). 



3 Random Matrix Theory for Supergravity 

In this section we will briefly review a few important concepts and results from random matrix 
theory, in order to make our analysis more self-contained. An accessible and fairly recent 
introduction can be found in [10]; see also the text by Mehta [TT] . 



3.1 Classical ensembles 

A foundational idea in random matrix theory is that given only limited statistical information 
about the entries of a diagonalizable N x N matrix, for large N one can make incisive state- 
ments about the statistical properties of the eigenvalues. For our purposes, the properties of 
principal interest are the eigenvalue spectrum for a typical matrix, and the probability of a 
large fluctuation of the smallest eigenvalue. 

We begin by reviewing the ensembles relevant for this work. 
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3.1.1 The Wigner ensemble 

One of the simplest and best-known ensembles of random matrices is the Wigner ensemble of 
Hermitian matrices, also referred to as the Gaussian Unitary Ensemble fTIj [T3| [H]. Elements 
of this ensemble, which we refer to as Wigner matrices, are N x N Hermitian matrices M given 
by 

M = A + A\ (3.1) 

where Aij ioi i,j = 1,. . . ,N are i.i.d. variables drawn from Q{0,a), and the dagger denotes 
Hermitian conjugation. 

The measure on the space of matrices is 

dP{M) = H f{M,,) dM,, , (3.2) 

l<i<jr<Ar 

where f{Mij) denotes the probability density of observing Mij. For normally-distributed entries 
of M, the joint probability density of the eigenvalues Ai, . . . , Aat is obtained by a unitary change 
of coordinates, 

N N 

/(Ai, . . . , A^) = C exp(-- 5^ A,2 + 2 ^ ln|A, - A,|) , (3.3) 

i=l i<j 

where C is an A^-dependent normalization constant. As conceived in the famous work of Dyson 
[T3] . this joint probability density can be given a physical interpretation in terms of a one- 
dimensional Coulomb gas of charged particles executing Brownian motion under the influ- 
ences of a confining quadratic potential and of mutual electrostatic repulsion. This physical 
picture has proved to be very fruitful in deriving exact results for a variety of properties of the 
eigenvalue spectrum (see e.g. [HE]), and in ^and ^we will see that repulsion between pairs 
of eigenvalues significantly impacts the stability of critical points in supergravity. 

At large A^, the eigenvalue spectrum of a Wigner matrix converges to the celebrated Wigner 
semicircle law, 

p(A) = ^y4iV^^. (3.4) 

where p is the eigenvalue density. Setting a = the eigenvalue spectrum has support in the 
interval [—2,2], cf. Figure [l} 



3.1.2 The Wishart ensemble 

The second class of random matrices we will need are complex Wishart matrices, which take 
the form 

M = AA\ (3.5) 

where A is an N x Q complex matrix with entries drawn from Q{0, a), and Q > N. The study 
of this ensemble dates back to Wishart's investigation of sample covariance matrices [16], and 
the universality evident in the Wishart ensemble provided some of the inspiration for Wigner's 
subsequent development of random matrix theory. 
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-2-1012 1234 

Figure 1: The eigenvalue spectra for the Wigner ensemble (left panel), and the Wishart ensem- 
ble with N = Q (right panel), from 10^ trials with N = 200. 



As a Wishart matrix is the Hermitian square of another matrix, it is necessarily positive 
semidefinite. The joint probability density of a complex Wishart matrix is (cf. e.g. (lOj) 



N N N 



/(Ai,...,Ajv) = C exp(-^^A, + 2 ln|A, - A,| + (Q - iV) ^InA,) . (3.6) 

i=l i<j i 

In the Coulomb gas picture, the non-negativity of a Wishart matrix corresponds to the presence 
of a hard wall at A = 0. 

The eigenvalue distribution in the Wishart ensemble is given by the Marcenko-Pastur law 
[T8] . which takes the form 

for the special case N = Q that will be relevant in our analysis, cf. Figure [1} 

The probability density function of the smallest eigenvalue Ai was first computed by Edel- 
man [T7], and for our purposes it suffices to note that for = Q and cr = its average 
position (Ai) scales as 

3.1.3 The Altland-Zirnbauer CI ensemble 



The matrix Ai appearing in the critical point equation (2.4) has an eigenvalue spectrum that 
is broadly reminiscent of the Wigner semicircle law, but the 2A^ eigenvalues of Ai come in 
opposite-sign pairs ±Aa, with < Ai < . . . < Aat. As observed in [2], matrices Ai of the form 



(2.5) belong to the Altland-Zirnbauer CI ensemble [19]. For normally-distributed entries of Ai, 



the joint probability density of the eigenvalues is 

N N N 



/(Ai, . . . , A^) = C exp(-^ ^ A^ + ^ ln|A2 - A^] + ^ In |A,|) . (3.8) 



a 

1=1 i-^j 
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-2 -1 1 2 - 0.()()5 0.000 0.005 0.010 

Figure 2: The eigenvalue spectrum for the CI ensemble, from 10^ trials with = 200. The 
full spectrum appears in the left panel, while the right panel shows the details of the cleft at 
A = 0. Notice that the boundary of the linear regime occurs for A ~ 4. 



In the Coulomb gas picture, the additional term J2iLi can be interpreted as encoding a 
repulsive force between each mirror pair of eigenvalues, ±Aj |I9]. This repulsion is particularly 
important for the smallest eigenvalue of A4, and leads to a linear cleft in the eigenvalue spectrum 
for small A: 

p(A) ^ A;A + C(A^) , (3.9) 

with k a constant of order unity, so that the eigenvalue density vanishes at A = 0, cf. Figure [2] 
Recalling that the critical point equation (2.4) requires that has 2\W\ as an eigenvalue, we 
see that critical points with very small \W\ are rare in comparison to those with 2\W\ ~ 1. 



3.2 Fluctuations of extreme eigenvalues 

The eigenvalue spectra presented above describe the typical configurations of eigenvalues: for 
example, the eigenvalue spectrum for an ensemble of Wigner matrices with entries drawn from 
f2(0, a) is zero outside [—2\/N ct, 2\/N a] , but this does not mean that no matrix in the ensemble 
has eigenvalues outside this range. Instead, the Wigner spectrum has a 'soft edge' at each end of 
the semicircle: there is a nonzero probability that one or more eigenvalues can be found beyond 
this edge. In contrast, the Wishart spectrum has a 'hard edge' at A = 0, as the matrices in 
question are necessarily positive semidefinite, while the other edge of the spectrum is soft. 

As we shall soon establish, the spectrum of the Hessian matrix "H extends to negative 
values, so that the presence of tachyons is generic, but not guaranteed. It is therefore essential 
to determine the probability that the smallest eigenvalue of "H happens to be large enough so 
that Ti is positive definite. We refer to this occurrence as a fluctuation to positivity. 
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3.2.1 Eigenvalue fluctuations and the Tracy- Widom law 

The study of fluctuations of the smallest (or largest) eigenvalue was initiated in the pioneering 
work of Tracy and Widom [3]; see also [20]. The theory of fluctuations is best-developed 
when the fluctuations are suitably small, with the deviation of the extreme eigenvalue from its 
mean position being 0{N~^/^) for the case of the Wigner ensemble. In this case, the smallest 
eigenvalue Ai is given by (cf. [21] for a useful summary) 

+ N-^'\, (3.10) 

where for — )■ oo, x is a random variable that follows the Tracy- Widom distribution F2 [3]. 
Extension of the Tracy- Widom law to the Wishart ensemble was achieved for real matrices in 
[22] . and for complex Wishart matrices in [23]. 

In ^we will find that for a typical supergravity critical point, a fluctuation of size 0{N^^/^) 
of the smallest eigenvalue of % is insufficient to render % positive definite. We therefore 
require an extension of the Tracy- Widom theory describing large fiuctuations, with the deviation 
from the mean position being as large as 0{-\fN). The theory of large fiuctuations has been 
developecj^ in a series of works by Majumdar and collaborators [H |26l [5l |271 [281 121], which we 
now briefiy review, focusing on the Wigner ensemble. 

Through a saddle point computation of the partition function in the Coulomb gas model. 
Dean and Majumdar [H |5] were able to evaluate the probability of a large fiuctuation of the 
smallest eigenvalue Ai of a Wigner matrix to the right of its mean position (Ai) = —2\/N. The 
result, at leading order in large A^, is [H [5] 



p(Ai > oc exp -2N^ip_ (t/Vm^ , (3.11) 



for t > —2yN and t + 2yN ~ 0{yN). Although ip-{y) is known in closed form, we present 
here only the result relevant for a fiuctuation to positivity: 

^*-(0)='f^. (3,12) 

(The corresponding result of [U [5] for real Wigner matrices agrees very well with the earlier nu- 
merical results of Aazami and Easther [6] for the same ensemble.) In summary, the probability 
that an A^ X A^ complex Wigner matrix is positive-definite is given by [HI HI |5] 



P oc exp 



(3.13) 



with c ~ At large A^ this is exceptionally small compared to the estimate P ~ 2^^ 

that follows from the naive assumption that the A^ eigenvalues are independent. Of course, 
eigenvalue interactions are fundamental to random matrix theory, so it is no surprise that 
omitting these interactions gives an entirely inaccurate result for the probability of positivity. 



^Earlier work on related fluctuations appears in |24j . For applications to counting critical points of random 
functions, see e.g. j25) . 
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An intuition from the Coulomb gas model will be helpful in our analysis. Consider a 
fluctuation of the smallest eigenvalue Ai to roughly the midpoint of the spectrum, as would 
be required for a fluctuation to positivity in the Wigner ensemble. The distance involved 
is 0{\/N), and 0{N) eigenvalues need to be displaced. As these eigenvalues experience a 
quadratic potential, the total energetic cost is O^N"^), consistent with the detailed results of 
in [5] . Similar results have been obtained for inward fluctuations of the soft and hard edges of 
the Wishart spectrum in [26] and [29] , respectively. The lesson is that a substantial inward shift 
of one or more eigenvalues has a statistical cost ~ exp(— A^^), and is hence extremely unlikely 
at large A^. 

Having assembled the necessary tools, we now turn to studying the stability of the Hessian. 



4 Stability of Generic Critical Points 



In this section we will study the Hessian "H at a generic critical point, where 



u, 



abc ■ 



(4.1) 



We refer the reader to Appendix |X] for a detailed demonstration that such points are indeed 
generic. 



In ^4.1 we examine the decomposition (2.12) of the Hessian matrix into a sum of terms 



and argue that each term is well-approximated as a member of one of the classical ensembles 



reviewed in ^ 3.1 We then obtain the eigenvalue spectrum analytically from the free convolution 



of the constituent spectra. In ^4^ we argue that the probability that a given critical point is 
a metastable vacuum can be obtained by adapting the results of [U |5] to the free convolution 
model. We then perform an extensive numerical analysis of the full Hessian matrix, finding 
that a generic critical point is exponentially unlikely to be a metastable vacuum. Thus, despite 
the abundance of critical points, this region of the random supergravity landscape is indeed a 
wasteland. 



4.1 The Hessian spectrum from a free convolution 



The Hessian T-L can be decomposed according to ( 2.12 ) as H = Hsusy + "^purc + T^k^^) + T^K^'-i) + 

"Hshift • 

From ^we immediately recognize that each of "Hsusy and ^^^{s) is very similar to a double 
copy of an A^-dimensional complex Wishart matrix. The correspondence is imperfect because 
Z (respectively K^^'>) is symmetric, so that the diagonal entries have twice the variance of the 
off-diagonal entries. We have verified that this minor difference does not significantly affect the 
eigenvalue spectrum. 

Next, Hpure and 7/^(4) can be modeled as 2A^-dimensional Wigner matrices. Once again, 
this is an approximation: the actual Hessian matrix at a critical point must incorporate the 



critical point constraint (2.4). The sum of two Wigner matrices is again a Wigner matrix, so 
we may write "Hpurc + ~ Wigner. 
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Assembling the pieces, and noting that the effect of 7/shift on the bulk of the eigenvalue 
spectrum is simply a translation, our model amounts to 

n ^ Wigner('Hpurc + ^k(4)) + Wishart('Hsusy) + Wishart(?^K(3)) • (4-2) 
4.1.1 Free convolutions and sums of random matrices 

To obtain the spectrum, we need to address a problem of the general form: "if A and B are 
random matrices with known eigenvalue spectra /i^, /^b, what is the spectrum [ia+b of their 
sum A+B?" If A and B were to commute, Ha+b would be the convolution of ha and /^s, but 
there is no justification for this assumption in our case. The solution of the general problem is 
provided by Voiculescu's theory of free probability. We will describe here only the immediately 
relevant tools of free probability, referring the reader to the text p] for details and references. 
Given two ensembles A, B of random matrices, the free convolution ffl is defined such that 

yUA ffl /is = ^^A+B , (4.3) 

i.e. the free convolution of the spectra of the summands is the spectrum of the sum. Just as 
cumulants are additive under ordinary convolution of random variables, free cumulants can be 
defined with the same additivity property under the free convolution. In principle Ha+b can 
be obtained from the R-transform, which is the generating function of the free cumulants [50] . 
but inversion of the R-transform can be rather cumbersome. For the large class of algebraic 
random matrices [31], which includes sums of Wigner and Wishart matrices, a more efficient 
approach [181 EI] is to work with the Stieltjes transform [31] . 

The Stieltjes transform of a probability measure d/i(a;) = p{x)dx with support on a real 
interval / is defined by 

,„,,(.) = / , ,4,4) 

Ji z-x 

where lm.{z) > 0. Algebraic random matrices are random matrices for which m^{z) is the 
solution of a polynomial equation in and z, e.g. the Stieltjes transform of the Wigner 
density solves the equation 

+ a zm^ + a = , (4.5) 

where a = (A^o"^) ^ ■ The probability density is readily obtained from m^j_{z) using the Stieltjes- 
Perron inversion formula, 

p{x) = — lim Im mJx + ie) . (4-6) 

Edelman and Rao have shown that the free convolution can be implemented efficiently through 
manipulations of polynomials involving the Stieltjes transform [31] . 



4.1.2 The spectrum as Wigner ffl Wishart ffl Wishart 



In terms of the free convolution ffl defined in equation (4.3), we may write the eigenvalue 
spectrum piTi) as 

p{n) ^ p(Wigner) ffl p(Wishart) ffl p(Wishart) , (4.7) 
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Figure 3: The histogram shows the spectrum of eigenvalues of the full Hessian matrix "H (2.12) 
for = 200 and u = .1, in units of F^, while the curve gives the analytic result (4.8) from the 
Wigner ffl Wishart ffl Wishart model, with no adjustable parameters. 



where p(Wigner) is given in equation ( 3.4[ ) and p(Wishart) is given in equation (3.7). Obtaining 
the Stieltjes transforms of p(Wigner) and p(Wishart) and using the polynomial method of 
we find the spectrum 



3u;4 + 117 + 9 {io^ - 4) A + 9A2 - (|)'/' ^{Xf/^ 



P(A) 



22/3 35/6 ^ ^(A)l/3 (c^2 + Q) 



(4.8) 



where cu 



V3\W\ 



F ' 



ip{X) is given by 



9uj^{X + 1) + 27a;2 (A^ - 5A + 24) - + 18A''^ - lOSX"^ + 216A + 1314 , (4.9) 



and 



T = -(^2 + 6)^ 81A^ + 162A^(tj2-7) +9A^(l3a;'^-180a;2 + 62l) + (4.10) 

+ 18A {2cu^ - 35co^ + 333^^ - 630) + Acu^ - 48a;^ + 873a;^ - 126366^^ - 9396 . 

This is one of our primary results. Figure [3] illustrates the remarkably good agreement between 
(4.8) and simulations of the full H. 



4.2 Eigenvalue fluctuations and de Sitter vacua 

Although we now have an analytic result for the eigenvalue spectrum in the Wigner ffl Wishart 
ffl Wishart model, which gives an excellent approximation to the spectrum of Ti itself, com- 
puting the probability of a large fluctuation for this model is rather involved. The literature 
summarized in §3.2| contains detailed characterizations of the fluctuations of extreme eigenval- 
ues of Wigner or Wishart matrices, but a direct computation of the large fluctuations in the 
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Wigner ffl Wishart ffl Wishart model would require a dedicated saddle point analysis along the 
lines of in |5], and is beyond the scope of the present work. 

From the Coulomb gas model, it is clear that a sufficiently large fluctuation of the smallest 
eigenvalue will be sensitive to the global shape of the spectrum: such a fluctuation will displace 
the eigenvalues to its right, with an energy cost that depends on their density. However, a small 
fluctuation will displace only the eigenvalues near the edge of the spectrum, and the likelihood 
of such a fluctuation should therefore depend only on the edge shape. Correspondingly, it has 
been conjectured pij that Tracy- Widom fluctuations will be seen in essentially any algebraic 
random matrix whose eigenvalue density has square root behavior at its edge. As our Wigner 
ffl Wishart ffl Wishart model falls in this class, we expect that small fluctuations of the smallest 
eigenvalue of H will be governed by the Tracy- Widom law. 

Examining the position of the left edge in the Wigner ffl Wishart ffl Wishart model, we see 
that a fluctuation to positivity is not a small fluctuation in the sense of [3] : the distance to the 
origin i^(9(A/iV), not 0{N~^^^). However, between the left edge and the origin, the spectrum 
has a shape very reminiscent of the semicircle law. Emboldened by this, we anticipate that the 
probability of a large fluctuation of the smallest eigenvalue of Ti is accurately modeled using 



the corresponding probability (3.13) in the Wigner ensemble, i.e. we expect P oc exp(— cA^^) 
with p ~ 2. 

4.3 Numerical results 

We now report on the results of extensive simulations of fluctuations to positivity in the full 
T-L model. These simulations make no approximation. We include the full structure of H 
(e.g., the slight difference between Hsusy and a Wishart matrix), and we do not rely on any of 
the analytical results reviewed in ^ No expansion in or in F is used. We simply create 
an ensemble of realizations of "H, following the prescription of ^ and directly determine the 
fraction that are positive deflnite. 

Naturally, these simulations could still fail to yield an accurate picture of the positivity 
probability in the supergravities derived from string theory: in particular, our deflnition of a 
random supergravity could be non-representative. Moreover, computational cost imposes an 
upper limit on N, and our extrapolation to larger N could be inaccurate. 



One detail of Ti requires further explanation. The critical point condition, equation (2.4), 
enforces that A4 has eigenvalue 2|iy|. For a given W, no matrix drawn randomly will have 
a precisely correct eigenvalue (this reflects the fact that critical points are a measure zero 
subset of all points). To impose this constraint, we note that if Ai has an eigenvalue Avi^ G 
[(2 - 2e)|Vr|, (2 + 2e)|ll^|], but Xw 7^ 2|W^|, the distortions of the spectrum compared to that 
found at a genuine critical point will be of order e. In our simulations, we have taken e = 10~^. 
In gthis issue will pose a greater challenge: for \W\ ^ 1/A^, a very small fraction of randomly 



drawn Ai will fall in [(2— 2e)|Vr|, (2+2e)|Vr|], cf. equation ( 3.9[ ), and it becomes computationally 



^For ease of comparison to Sj, we take il = 7V(0, 1) in this discussion, although we set fl = Af{0, ^) 
elsewhere. 
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costly to find examples. 

Figure |4] presents the results of simulations of the full Hessian matrix 1-L (upper curve) and of 
the Wigner ffl Wishart ffl Wishart model (lower curve). The qualitative properties are similar, 
but the best-fit values of p are somewhat different. This is not surprising, as the numerically- 
accessible values of are not large: for > 20, stability is extremely rare, and it is difficult 
to obtain sufficient statistics to characterize the probability of positivity. For < 20, the large 

expansion underpinning our random matrix theory approach is marginal at best, with two 
important consequences. First, the correspondence between the model and simulation spectra 
is imperfect for A^ ~ 20 (contrast the superb agreement for A^ = 200 shown in Figure |3|, and 
the corresponding difference between the left edges of these spectra contributes to a different 
fluctuation probability. A primary cause of the difference between the spectra of the full Hessian 
matrix l-i and of the analytical model (4.8) is that T^pure involves the matrix Z, as does Hsusy, 
so that in the Wigner ffl Wishart ffl Wishart model, the Wigner matrix is correlated with one of 
the Wishart matrices. For A^ ^ 1 (and also for u <^1) this correlation becomes less important. 
Second, for small enough A^, fluctuations to positivity are governed by the Tracy- Widom law 
(3.10) rather than by the considerably steeper large-fluctuation expression (3.11), so that a fit 
that includes data points starting from N = 2 will result in a value of p that is smaller than 
the asymptotic large A^ value. 

In summary, due to the challenges inherent in studying extremely rare events numerically, 
we have not obtained sufficient data at large A^ to make a definitive determination of the large 
A^ behavior of the probability, and this is an interesting problem for the future. In light of the 
arguments of §4.2 it remains reasonable to conjecture that p ~ 2 at sufficiently large A^. 

Next, Figure [5] shows the trends in c and p, cf. equation (1.1), as to is varied in [0.1,1]. 
There is a distinct increase in c as w increases, while p shows a barely significant increase, so 
that overall the probability of positivity drops substantially as u increases. This trend can 
be understood as follows: increasing u at fixed F reduces Hshift, and hence shifts the entire 
spectrum toward more negative values, making a fluctuation to positivity more improbable. 

We note that for fixed F, increasing u reduces the cosmological constant, so within this class 
of critical points, the probability of stability is higher at higher cosmological constant. One 
should not read too much into this, however, as one can increase the cosmological constant by 
increasing F, \ W\, msusy by a common factor without affecting the probability of stability. 



5 Stability of Critical Points with Small F-terms 

The conclusions of ^ apply to a generic critical point, by which we mean one at which the 
functions K, W are random functions that do not automatically manifest any special hierar- 
chies!^ However, it is far from clear a priori that a typical metastable vacuum arises from 
among the set of generic critical points: a tiny subclass of critical points that enjoy a high 

^Of course, large ratios can arise in this setting by chance, but this possibility is already encoded in the 
results of ^ 
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Figure 4: The logarithm of the probabihty P(Amin > 0) that the smallest eigenvalue of "H is 
positive, as a function of N, with u = 1. Upper branch: simulations of the full Hessian matrix 
T-L, with best-fit values p = 1.50 ± 0.10, c = 0.29 ± 0.06. Lower branch: simulations of the 
Wigner ffl Wishart ffl Wishart model, with best-fit values p = 1.90 ± 0.04, c = 0.21 ± 0.02. The 
error bars give the 2a statistical uncertainty. 




0.10 



0.15 



0.20 



0.25 



0.30 



0.35 



Figure 5: The ellipses show the 2a allowed regions of the p — c plane, cf. equation (1.1), for 
u = 0.1, 0.2, . . . , 1.0, from left to right, with 2 < < 23. As u increases (so that for fixed F 
the cosmological constant decreases), c increases substantially, while p increases slightly. 
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likelihood of stability as a consequence of some special structure might well account for most 
of the metastable vacua. 

As originally noted by Denef and Douglas in [2], a particularly interesting class of critical 
points are those at which the F-terms are small compared to the supersymmetric masses: 
approximate supersymmetry can be expected to make stability more likely. Specifically, we 
will consider de Sitter critical points at which 

Vl\W\ < F < \Zab\ ~ \Uabc\ ■ (5.1) 

In this section we will reexamine the stability of critical points in this approximately-supersymmetric 
corner of the supergravity landscape. As will become clear, our conclusion differs from that of 
Denef and Douglas, and we will carefully explain the reason for the disparity. 

We will see that eigenvalue repulsion in the mass matrix between the bulk of the eigenval- 
ues and the Goldstino direction typically generates at least two tachyonic directions, rendering 
generic critical points unstable. This effect significantly influences the fine-tuning needed to 
obtain a metastable de Sitter solution in supergravity. Through numerical simulations and 
through statistical analysis we will show that metastable critical points constitute an exponen- 
tially small fraction of all critical points: A^vacua — A'crit.pts.e"'^^*', with p ^ 1.3 and c ~ 0.08. 



5.1 The Denef-Douglas landscape of de Sitter vacua 

In order to analyze the stability properties of the mass matrix in the regime where F <^ mgusy, 
we write "H as [2], 

U = v;; + V^ + V^, (5.2) 



with 



K'= ( ° ) ' S, = UabcF\ (5.4) 



and 



V^' = {M + \W\1){M-2\W\1) , (5.3) 

Si 
Si 



= ( I ) , S, = 6,iF' - F,F-, - Ki^^F^F'^ , (5.5) 



where M. is the matrix given in equation (2.5) 



When msusy ^ F, all but two of the eigenvalues of 1-L are generically of order ml^^^, and are 
predominantly determined by V^", with corrections of order F from V^', and of order from 



V2 ■ However, the critical point condition, equation (2.5), requires that M. has an eigenvalue 
2|VF|, with the corresponding eigenvector pointing in the Goldstino direction. The eigenvalues 
of M. come in pairs differing only in sign, and the eigenvalues \± = i:2\W\ of Ai correspond to 
eigenvalues ml = 0, ml = 4|Vrp of Vq. The larger of these 'Goldstino' eigenvalues is C(F^), so 
that one cannot a priori neglect the effects of V{', on the stability of the Goldstino direction. 
This section is dedicated to a careful examination of these effects. 
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Setting Fa = d^Fe^'^'' and performing a unitary transformation that diagonalizes ZaZi^^ 
we obtain the simphfied mass matrix 







ab 


Kb 






ab 


Kb 



( 



0{F') diag(A^,) 


Kb 


Kb 


0{F') diag(A^,) 



(5.6) 



The eigenvalues of Z^Zi^ have been denoted A^,, for a' = 2, . . . A^, while by the critical point 
equation we have m^j = 2|iyp — i^injF^. In the approximately-supersymmetric regime, A^/ ^ 
TTisusyF, and we have correspondingly omitted 0{F'^) contributions to the diagonal entries 



{Hi,) ^1^1 for a' = 2, . . . N. The notation emphasizes that the matrix appearing in (5.6) is 
not a truncation of H to some order in F. Instead, has been strategically simplified so that, 
while it efficiently yields results for the two smallest eigenvalues of H that are accurate up to 
0{F^) corrections, the higher eigenvalues of "H*, which are generically positive in any case, do 
not coincide with those of "H to this accuracy. 

Following the discussion of [2] , we focus on the submatrix spanned by the normalized eigen- 
vectors of Ai with eigenvalues ±2|Vr|, which in the above basis can be expressed as 

= ^(e^^'^ 5/ + e-^^X^-^') and (M^rO, = ^{e'^' - e^^'' 6f^') , (5.7) 

where A-d = -dp — "^w- Neglecting for the moment mixings with the other eigenvectors of Ai, 
this 2x2 Goldstino submatrix of the full mass matrix is 



n 



sub 



m 
m 



m 



11 



2_ 

11 

II Ki 



(5.8) 



While a diagonalization of this subsystem by itself does not in general correspond to a diago- 
nalization of the corresponding directions in the full mass matrix, it is instructive to attempt 
to treat the off-diagonal mixings in perturbation theory. The eigenvalues of the submatrix are 
given by 

= mlj ± \ml^\ , (5.9) 



which can be written as 

= 2\W\'^ - RniiF'^ ± 



UiiiFe-'^ 
l-Riiiil 



2\W\'e 



(5.10) 



0(F°) is the term UmF, so that 



The dominant contribution to for \Uui\ 
generically h~ < 0. 

In [2], it was observed that fine-tuning |f/iii| to be 0{F) is necessary for stability of the 
mass matrix. However, [2] also argued that this condition is sufficient, and concluded that 
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metastable critical points are fairly common in a supergravity landscape. We will now show 
that the eigenvalues /i^ of the submatrix Tisuh of Ti cannot be regarded as good approximations 
to the actual eigenvalues of Ti. Upon computing the leading-order corrections to (5.9), we will 
find that metastable critical points constitute an exponentially small fraction of all de Sitter 
critical points. 



5.2 Eigenvalue repulsion induces tachyons 

In this section we discuss the correction to the eigenvalues of the mass matrix induced by V^' , 
equation (5.4). To determine the two smallest eigenvalues to (9(-F^), we may neglect (9(-F^) 
contributions to fn'^,-j^, for a' ^ b'. With this simplification, equation (5.6) can be written as 



Hi, 



1 







S 


T 


\ 







diag(A2,) 


Va' 


Ta'b' 






s* 




















diag(A2,) 


/ 



(5.11) 



In equation (5.11 ) we have introduced the U{N — 1) scalar s 
the symmetric tensor Ta'b' = "^a'b'- The unitary transformation 



m 



11) 



the vector Va' = ml^,, and 



U 
















^a'b' 








1 





1 














^a'b' 



with a = arg(s), diagonalizes the sub-matrix T^sub, i-e. 



(5.12) 



( ^' 


ui' 





T 


\ 


Ua' 


diag(A2,) 


Ua' 


Ta'b' 







4 




-Wa' 






^a'b' 


* 

-Wa' 


diag(A2,) 


/ 



(5.13) 



where u„i 



-^J^Ua', Wa' 



a', and /i^ is given by equation (5.9). The conclusion of [2] 
is that modest fine-tuning of s ensures the positivity of h^, and hence of T-L. Here we investigate 
the effect of the vector Va' on the eigenvalues of Ti. The leading-order effect of the tensor Ta'b' 
is to induce 0{msusyF) shifts of the eigenvalues A^, ^ mgusy-^, so that we may consistently 
neglect Ta'b'- The characteristic polynomial of is then given by 



C(P) 



+ 



Co(p) 

N 

E 

a',b'=2 



1 - 



N 

E 

b'=2 



\Vb' 



N 



{h^ - p){Xl, - p) 



Va' Vb' 



-E 

b'=2 



\Vb' 



{h- - p){\l, - p) 



{h+ - p){h~ - p){\l - p){\l - p) 



(5.14) 
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where Co(p) denotes the characteristic polynomial for Va' = 0, i.e. 



N 



Co(p) = {h- - p){h^ - p) n i^l' -pY- (5-15) 

a'=2 

The leading-order effect of the vector Va' is evidently to induce an 'interaction' between the 
approximate eigenvalues and A^,. This interaction is a manifestation of eigenvalue repulsion, 
and indeed, the effect of each term in the sum is to increase the splitting between and X^,. 
Restricting the polynomial to small values of p close to the smallest eigenvalues of the mass 
matrix and dividing by the overall factors of the larger eigenvalues, equation (5.14) can be 
rewritten as 

^/ X ^ \ \2 ^ \ \2 

- (A* - P){h- - P) - (h- - E ^ - - E ¥ 



nAT ,0 ^ A^/V" V" \2 y ^' \2 

a'=2^a' b'=2 b' ^,'=2 b' 



N 

+ E ("6) 



\Va'r_\Vb' 
a', b'=2 ^^'^b' 



The smallest eigenvalues of the mass matrix are thus given by 



Af I 19 N 

9 , -I- ' 

m 



2 12 



b'=2 ^b' fe,=2 ''b' 

b'=2 ^b' 

This is one of our main results. The smallest eigenvalue of the Hessian matrix Ti differs 
from that of [2j by the non-positive term —-^^^^=2 ^^la''^ ; in a manifestation of eigenvalue 

b' 

repulsion between the Goldstino and the supersymmetrically stabihzed moduli with masses of 
order A^,. 

We now turn to assessing the impact of this contribution to m^. 



5.3 Eigenvalue fluctuations and de Sitter vacua 

In this section we will determine the probability that a randomly chosen approximately- 
supersymmetric critical point is metastable by computing the probability that m^, as given in 
equation (5.17), is positive. Our approach is to determine the statistical propertied of each 
term in the sum, i.e. we will obtain the cumulative distribution function (cdf) for each term, 
from which the corresponding probability density function (pdf) can be obtained by differen- 
tiation. Although in principle one might hope to convolve the constituent probability density 



^Recall from [2] that a fine-tuning of |C/iiii^| < is necessary for stabifity, and granting this fine-tuning, 
all the terms in equation (5.17) are of the same order, 0{F^). 
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functions to obtain the pdf of m^, this is rather involved. Fortunately, we will find that one of 



the terms of equation (5.17) dominates both in magnitude and in the probability of fluctuations, 



and it suffices to examine this term. 

For the analytical estimates provided here, we will assume fl = A/'(0, ^), though similar 
arguments could be made for e.g. the uniform distribution. 



We find it convenient to rewrite (5.17) as 



m 



with 



and 



where 



r=-uj^ + K,,-Kne - Kun - \t 



hoi 



s 



N 

E 



l^llb'P 



ol 



U,,,e^^^^-^^^^F~' - -co' 



and we have used the definition (2.30). 



(5.18) 
(5.19) 

(5.20) 

(5.21) 



5.3.1 Subdominant contributions 

We will begin by studying the terms collected in T, which do not involve the vector Va'- 

At the critical points of interest, \/3|l^| < F, so that the total energy density is nonnegative. 



Thus, the first term in equation (5.19) gives a contribution in the range [0, 3]. 

The second term of equation (5.19) is jK'^'^^p = K^^^Kn^, which is the sum of squares of N 
random variables, each drawn from M{0, -^). Thus, is distributed as -^Xn^ where Xn 

a chi-square distribution with N degrees of freedom. Since Xn mean N, we conclude that 



(|i^(3)|2) = l. 

To find the probability of fiuctuations, we note that the corresponding cdf is given by 



< X 



P{^X% < .) 



r(iv/2) 



7 



Nx 



2 / 



(5.22) 



(5.23) 



where 7 denotes the lower incomplete gamma function. The asymptotic behavior can be ob- 
tained as follows: by the central limit theorem, the cdf of a chi-square distributed variable for 
^ 1 degrees of freedom tends to that of a Gaussian distributed variable with unit variance, 
P{X% < 2/) ~ ^(^(0, l)<x), where x = ^ [32]. 

We are particularly interested in the probability of iJC'^^^p fiuctuating to a large value and 
thereby stabilizing the smallest eigenvalue oi T-L. As we will describe below, large in this 
context means 0{N), so that we consider 



P(|ir(3)|2 < ^ p(a/'(0, 1) < 2-1/2^3/2^ 



(5.24) 
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for ^ 1, from which we obtain 

P(|i^(3)| 



> N] < 



Ar3/2 



(5.25) 



Fluctuations of jii'^^-'p may therefore be neglected in comparison to the much more probable 



fluctuations we will discuss in ^5.3.2 



The third term in (5.19), -ft'ljl^j, is normally distributed with a vanishing expectation value, 



and with a variance no larger than so that large deviations of the order N are likewise so 
improbable as to be negligible: 



P{k[\\-^>N] ~e 



-Ar3 



(5.26) 



The fourth term in (5.19), — |thoi|; is negative semidefinite, and only one entry (not eigen- 



value) of namely f/m, needs to be adjusted in order to change the size of |thoi|- Therefore, 
it is straightforward to fine-tune |thoi| to be small. It is clear from the discussion above that, as 
originally noted in [2], m?_ is generically negative unless Um is fine-tuned to make Ithoil ^ C'(l). 
For our goal of obtaining a conservative estimate of the probability that > 0, it suffices to 
set Ithoil = 0. 



5.3.2 The eigenvalue repulsion term 



Finally, the last term in equation (5.18) is the sum of squares of — 1 terms. The numerators 



of the terms in equation (5.20) are the squares of independent normally distributed variables. 



while the denominators are the squares of the eigenvalues of Ai. 

The eigenvalues of range from around 0{j^) to 2 in unit^of msusy, so {\S\) ~ 0{N). 
Recalling that the contributions to T have mean sizes independent of A^, we conclude that S 
provides the dominant contribution to at large N. Since (S) < 0, this term destabilizes 
generic critical points in the approximately-supersymmetric regime. 

To determine the (small) probability that is nevertheless positive, we will now estimate 
the probability that S > —1, so that T + S can be positive. First, we recognize that in light of 
the discussion in ^ fluctuations that increase the denominators appearing in S, corresponding 
to inward fluctuations of the eigenvalues of a Wishart matrix, are extremely unlikely at large N 
|29] . Fluctuations of S toward smaller magnitude are principally determined by the fluctuations 
of the numerators. (We have explicitly verifled this in simulations.) This justifles simplifying 
the problem by flxing the factors of A^, to their mean values, (A^,), as determined by the bulk 



distribution given by equation (3.7). Henceforth we consider the sum 



N 



lb' 



b'=2 



(5.27) 



'Since by assumption Uatc ~ Zat, the dependence on the supersymmetric mass scale rrisusy cancels between 



the numerator and denominator in equation (5.20) 
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which is the weighted sum of — 1 variables that are all independently distributed as Xi- 
Weighted sums of x^-distributed variables (or, equivalently, sums of F-distributed variables 
with different scale parameters) occur frequently in statistics, and in particular in the theory 
of the distributions of quadratic forms. While we have not found a closed-form expression 
for the convolution of — 1 such terms, approximations for expressions like (5.27) have been 



developed. An approximation by Solomon and Stephens has been argued to be particularly 
accurate in the small-argument regime of interest [33] , but we will find that for our purposes it 
does not constitute a close approximation to the cumulative probability for small arguments. 
By matching the first three algebraic moments /xi, fi2, A^s of S' to those of a ■ w^, where w is 
X^-distributed and a, 6, and r are constants, this approximation is obtainec^ by numerically 
solving the equations 



1^1 



a 2' 



r(fe + i) 
r(i) ' 

r\ r(26 + 



[Tib + 

2 r(3&+ 1; 
[r(&+i 



3 • 



(5.28) 
(5.29) 

(5.30) 



Thus, in this approximation, 

P{S' <s)^ P[a{xlf <s)= P[xl < ( 



a 



(5.31) 



As T ~ C'(l), 'irn?- could be positive if S' fluctuates down to be C(l), for which we obtain 



P mi > 



p[s'<i 



-c-NP 



(5.32) 



where c ~ 23, and p ~ 0.24. As we will see in §5.4 even though this approximation qualitatively 
matches the shape of S', for A^ ^ 1 it severely overestimates the probability of a fluctuation of 
S' to be of 0{1), and it remains an open question to obtain a good analytic or semi-analytic 



approximation of equation (5.27). 



5.4 Numerical results 

Figure [6] shows a histogram of and its constituent terms T and 5, for A^ = 40. It is 
clear that S gives the dominant contribution to m^. Moreover, the narrow support of the T 
histogram illustrates the finding of ^5.3.1 that large fluctuations of T are much less probable 
than correspondingly large fluctuations of S. 

Finally, Figure [7] presents the result of simulations of the mass matrix in the approximately- 
supersymmetric regime. (The value of u has a negligible effect on stability in this regime.) The 
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There is a misprint in the fifth equation of §3.1 of 
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Figure 6: Histograms of the smallest eigenvalue m^, cf. equation (5.18), and its constituent 
terms T and S, for = 40, in units of F"^. The eigenvalue repulsion sum S has the leftmost 
peak, the total mass has the central peak, and T appears on the right. Note that S, and 
consequently m^, has support over a range of size N (not fully shown in the figure), while T 
has variance j^. 



data agrees well with (1.1), with p = 1.28 ±0.03 and c = 0.083 ± 0.008]^ This is a much larger 
value for p than that obtained by the analytical estimate of [: 5.3.2 so that the latter gives an 
extremely conservative upper bound on the asymptotic large probability of positivity. 



6 Beyond Random Supergravity 

In this section we will discuss potential extensions of our assumptions ( ^6.1 ), explain the con- 



sequences of decoupling for the probability of positivity (^6.2), and illustrate our results in the 



example of the KKLT scenario (^6.3) 



6.1 Universality 



The results of ^and ^were obtained from the assumptions enumerated in §2.3[ most notably, 
K and W were taken to be random functions of N scalar fields, so that their various (appro- 
priately covariant) derivatives are i.i.d. variables drawn from a distribution a). Equivalent 
assumptions are standard in the statistical study of flux compactifications, cf . [T] , and in partic- 
ular are fully consistent with the assumptions of |2] . Nevertheless, in this section we will venture 
a few remarks about possible extensions of this simplest definition of a random supergravity. 

First, we have taken the random variables to be normally distributed, Q = A/'(0, ■^), 
throughout this work, and now we justify this assumption. The celebrated phenomenon of 



"^To obtain a conservative bound, we fit to the data points with N > 7. 
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In P(An,i„>0) - In P(A^i„>0 | N=2) 

— , — , — , — , — : — , — , — , — , — : — , — , — , — , — ^ N 
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Figure 7: The logarithm of the probabihty P{Xmm > 0) that the smallest eigenvalue of "H is 
positive, as a function of A^. Each point corresponds to 10^ realizations of the full mass matrix, 
and the error bars give the 2a statistical uncertainty. The curve shows the best fit to equation 
(1.1), with p = 1.28 ± 0.03 and c = 0.083 ± 0.008. The dashed line with p = 1 is for reference. 



universality in random matrix theory ensures that at large the eigenvalue spectrum, and 
also the fluctuations of extreme eigenvalues, are independent of the statistical details of the 
inputs. Universality has been demonstrated in a staggering array of physical and mathematical 
systems, many with < C(IO^), including interfaces in hquid crystals [M], the timing of buses 
in Cuernavaca [35] , and the power output of coupled lasers j36] . See [371 EH] for overviews of 
universality and [391 HQ] for results with close connections to ensembles studied here; extensions 
to ensembles in which the matrix entries have power-law tails include [H]. The lesson is that 
the particular choice of Q is immaterial, provided that the moments of fl are appropriately 
bounded. (One should compare distributions ^1,^2 yielding the same root-mean-square size 
for entries in "H, as this sets the physical scale.) 

Despite the strong expectation that universality should be applicable for our system, it 
is still reasonable to ask whether the values of N in our analysis are large enough for these 
asymptotic results to apply in practice. We have addressed this point directly by repeating our 
simulations for different choices of distribution, with excellent agreement. 

A more fundamental question is whether in the effective theories derived from string com- 
pactifications, the derivatives of K and W are accurately modeled as i.i.d. variables drawn 



from any distribution, or if instead these quantities are not i.i.d, Microphysical constraints, 



for example the relics of extended supersymmetryj^ might be expected to introduce correla- 
tions among these variables, as in the special geometry relation of (2.23), cf. [21 [7], so that the 
derivatives of K and W are not all independent. A definitive answer to this question is beyond 
the scope of this work, but it is encouraging that in simpler cases such as the Wigner ensemble. 



^^We thank M. Douglas for instructive correspondence on this point. 
^■^See |42j for related work in maximal supergravity. 
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universality has been shown to apply to matrices with highly correlated entries [43j . 

It would also be interesting to understand the possible impact of global constraints on our 
considerations. We have studied the ensemble of critical points arising in a general supergravity 
theory, and we have taken the derivatives of K and W evaluated at each such point to be 
random functions. To understand the distribution of vacua within the moduli space, one should 
incorporate further structure. The (index) density of supersymmetric vacua is well known to 
be correlated with the curvature of the moduli space [44j|, while the global structure of the 
superpotential is better modeled as a random holomorphic section of a line bundle over the 
moduli space (see [S] for a definitive treatment of the density of supersymmetric vacua in 
this context). Extending the study of non-supersymmetric vacua to this level of detail is an 
interesting problem for the future. 

One might expect that constraints from Morse theory will require some small deviations 



from the purely statistical results obtained here^^ The random matrix ensembles we have 
described (as in [2] and earlier works) predict certain ratios between the numbers of saddle 
points of varying index, which are not automatically consistent with the Morse inequalities. 
The necessary adjustments can be accommodated without changing the number of minima 
relative to saddle points, and we find it plausible that any effect on the relative number of 
minima can be neglected. 

6.2 Decoupling improves stability 

A significant assumption in our analysis is that W and K are general random functions of N 
scalars. In physically well-motivated examples, there can of course be two or more sectors of 
fields with distinct mass scales. For instance, consideif^ a two-sector supergravity theory with 
= Nh + Nl fields, in which the heavy scalars 0^, a = 1, . . . Nh, receive large supersymmetric 
masses mn, and supersymmetry is dynamically broken in a decoupled system of lighter scalars 
0}^, i = 1, . . . Nl, at a much lower scale uil- 

Explicitly, such a model can be constructed from a superpotential and a Kahler potential 
that are additively separable. In a convenient Kahler gauge, one has 

K{4>h,4>H,4>L,4>l) = Kh{4>H,4>h) + KL{(f)L,(f)L) , 

W{<PhAl) = Wh{<Ph) + Wl{<Pl). (6.1) 



By assumption Zab ~ ran and Zij mi, while by (6.1), the cross-couplings in the supersym- 
metric mass matrix are small: Zai = KaWi ~ 0{F). Thus, at small F, T^susy separates into 
two distinct Wishart matrices. (If F is not small compared to m^, or if the separability of the 
superpotential is imperfect, then the off-diagonal masses in T^susy cannot be neglected.) 

A cautionary remark is necessary at this point. The masses-squared in a supersymmetrically- 



stabilized sector are not necessarily positive: setting F = in equation (2.12), the contribution 



^^We thank B. Czech for very helpful correspondence about these constraints. 
^^We are indebted to S. Kachru for emphasizing the importance of this example. 



28 



of ?/susy is nonnegative, but Hshift and "Hpure make tachyonic contributions that will be signifi- 
cant unless W <^ rrisusy Of course, the resulting masses do obey the Breitenlohner-Freedman 
(BF) bound, but this in itself does not guarantee that this sector will remain stable after 
supersymmetry is broken badly enough to make the cosmological constant positive. 



that for F <^ rrisusy, superpotential couplings of the form 

A' 



We now recall from ^5 

contribute to the destabifization of the Goldstino direction, cf. equation (5.17), where A' runs 
over all fields. The numerator of this contribution to the Goldstino mass from the heavy, 
supersymmetric subsystem is 

\Uua\' = \VaZn\' = \daZn + K,Zn\\ (6.2) 



which under the decoupling assumptions of equation ( |6.1 ) is of order (For a non-decoupled 

system one finds instead |f/iiap ~ "^susv-) 



Since the denominator is of order 171%, the separability of equations (6.1) leads to a 
suppression of order of the heavy fields' negative contribution to the Goldstino direction 

rrijj 

mass-squared. Thus, even for a modest hierarchy between the supersymmetric masses, m'j^ > 
the high-scale sector decouples, and does not contribute significantly to the mass of 
the Goldstino. 

In conclusion, the relevant number of fields for the stability analysis of ^is Nl, the number 
of 'light' fields that participate in dynamical supersymmetry breaking (the superpartner of the 
Goldstino is assumed to be entirely among these fields.) The fraction of critical points that 
are metastable is then proportional to exp(— cA^^^). Provided that the heavy sector, taken in 
isolation, contains a number of supersymmetric vacua that is exponential in Nh, then the net 
result, for Nh S> Nl, is a mild reduction in the number of metastable vacua. 

6.3 Stability in the KKLT scenario 

The KKLT scenario |16] provides a useful setting to illustrate our findings. Consider a model 
with Nk = /i(^) Kahler moduli Tj, z = 1, . . . Nx, and Nc = h'^'^ complex structure moduli Ca, 
a = 1, . . . Nc- Suppose that the superpotential takes the form 

Nk 

I /TT \ 

(6.3) 



W= fcAn + Y,MC) exp (—T^ 



where rij is the dual Coxeter number for superpotential terms generated by gaugino conden- 
sation, rii = 1 for terms generated by Euchdean D3-branes, and j G A Q depends on the 
(a- Finding a compactification with many moduli for which each Kahler modulus appears in 
the nonperturbative superpotential is a difficult task (cf. [l?! HH] for detailed examples). Our 



purpose is to show that, granting a superpotential of the form (6.3), then for N^ ^ 1, an 
exponentially small fraction of de Sitter critical points are metastable vacua. 

An important scale in the problem is the flux scale maux, which sets the typical size of the 
supersymmetric masses for the (a- In light of the very large number of choices of quantized 
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flux [49j, one can find configurations in which the vacuum expectation value of the classical 
superpotential obeys {J G A Q) ^ "^1^^. This fine-tuning is necessary in order to obtain a 
parametrically controlled vacuum with a reasonably small cosmological constant. Given such a 
flux superpotential, one can find [IB] a supersymmetric AdS vacuum with all moduli stabilized. 

Our goal is to assess the stability of such a configuration after uplifting to de Sitter space. 
As a conservative first step, we imagine that the uplifting increases the cosmological constant 
without creating new instabilities, as a (fictitious) moduli-independent D-term would do. We 
expect that more plausible sources of positive energy will worsen any instability problems seen 
in this simple case. 

To begin, we will examine the masses in the supersymmetric AdS vacuum, and ask whether 



these masses can be made positive definite and large compared to |iy|^^ If they can, then 
after a rigid uplifting to de Sitter space (in the sense described above), the mass matrix will be 
dominated by Hsusy, which is positive definite. 

The dependence of Ai on the (a can be neglected self-consistently for Zab, but since the 
nonperturbative contributions lead to an imperfect separability of the superpotential, mass 
mixings through terms of the form Zai cannot be neglected, and the scale of the entries Zai and 
Zij is now \W\. Specifically, 

Zij = V,VjW = d,djW + {K,,-K,K,)W, (6.4) 
Zaj = VaVjW = dad.W - KjKaW , (6.5) 

where we have used the F-flatness conditions DiW = DaW = 0. The derivatives of K will 
not lead to enhancements in a controllable regime, while from (6.3) it follows that didjW ~ 
dadjW ~ W. 

The entries of "Hsusy with both indices in the complex structure directions are of order 
"^flux ^ while the mixed entries in "Hsusy receive contributions of order mfiux \W\ from 

terms of the form Z^Z-^^. The entries in the Kahler moduli directions are of order \W\^ . As in 



the discussion in §6.2[ the eigenvalues of "H split into two groups: the first consisting predom- 
inantly of the complex structure moduli and axiodilaton, which are stabilized at a high scale 
without BF-allowed tachyons, and the second consisting predominantly of the Kahler moduli, 
which have masses of order \W\. Since the supersymmetric Kahler moduli masses are not 
parametrically larger than the negative shift term "Hshift, or the off-diagonal contribution 7/pure, 
BF-allowed tachyons are typically abundant in the Kahler moduli sector in the supersymmetric 
AdS vacuum. 

Assuming a rigid uplifting to a de Sitter critical point, we recognize that the Kahler moduli 
sector constitutes a particular variant of the analysis of ^with F = 0: the somewhat more 
favorable regime described in ^ which requires ^ l^ijl) is inaccessible. Notice that in the 
generic regime of ^ the Goldstino is by no means the only tachyon, so that instabilities will arise 
in the Kahler moduli sector even if the Goldstino direction belongs to some other sector, e.g. a 
local dynamical supersymmetry breaking sector that engineers the positive vacuum energy. 
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We assume throughout that F < few x 
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We conclude that if a system described by a superpotential of the form (6.3) is rigidly 

uplifted to positive vacuum energy, the fraction of de Sitter critical points that are metastable 

1,1 
\+) 

'(+) 



vacua is proportional to exp{—cN^), with p > 1. For compactifications in which = /^('^ 
is not large, this is not a serious constraint, but it has significant impact for h^'\ ^ 1, and 



particularly for /i^^^ ^ /i^'^. 

One might object at this point that the arguments in favor of the existence of an approximately- 
supersymmetric regime as in ^ should hold for general functions W, K, so why are they 
not applicable here? The answer is simply that a superpotential of the form (6.3), which con- 
tains one single-instanton term for each Kahler modulus, is not a sufficiently general function. 
An obvious extension is to consider multiple terms (i.e., a racetrack) for each of the Tj. For the 
purposes of this discussion, we grant any topological prerequisites for such a multiple racetrack, 
e.g. we suppose that the compactification admits more than one stack of D7-branes in each ho- 
mology class. Then, by fine-tuning the fiuxes to adjust the prefactors Ai, cf. ^U\, one can 
plausibly arrange that the diagonal entries of Zij are large compared to \W\. However, there 
is a statistical price for this fine-tuning, of order {\W\/msusy)^'^ ■ Recalhng that the boundary 
between the regimes of ^ and ^occurs for |W^|/msusy ~ '^/Nk, this fine-tuning is of order 
A^-^^, which can be significant. 

To recap, if one assumes a rigid uplifting that changes the cosmological constant without 
changing the moduli mass matrix, then for a superpotential of the form (6.3), the Kahler 
moduli sector will have supersymmetric masses of order \W\, and will be governed by the 
instability analysis of ^ (with F = 0, W ^ 0), with positivity probability P oc exp(— c A^^), 
with p > 1. By fine-tuning a superpotential involving 0{Nk) racetracks, requiring a statistical 
price ~ A^^^^, one can make the supersymmetric masses large enough to guarantee stability. 

The situation is considerably worse if supersymmetry is spontaneously broken by an F-term 
in the Kahler moduli sector: again the analysis of ^ applies generically, but even after fine- 
tuning 0{Nk) racetracks as above, the Goldstino instability will still fall in the Kahler moduli 
sector, so that the instability analysis of ^ with p ^ 1.3, is apphcable. 

In summary, instabilities appear generic in the Kahler moduli sector after uplifting, with 
metastable de Sitter vacua constituting a fraction < exp(— c A^^) of all de Sitter critical points. 
However, for /i^'^ ^ ^{+)' number of KKLT vacua remains astronomically large, and the 
overall status of the model is not altered by our findings. 



7 Conclusions 

We have considered a general four-dimensional TV = 1 supergravity theory whose superpo- 
tential and Kahler potential are random functions of A^ ^ 1 scalar fields, and asked what 
fraction / of de Sitter critical points, with supersymmetry spontaneously broken by an F-term, 
are metastable vacua rather than unstable saddle points. Our conclusion is that an exponen- 
tially small fraction of critical points are vacua: / oc exp(— cA^'^), with p > 1.3, which differs 
significantly from earlier results implying f ^ jr- 
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The character of the instabihties that arise depends on the relative sizes of the supersym- 
metric and supersjTumetry-breaking masses. At a generic critical point, the soft masses are 
comparable to the supersymmetric masses, and supersymmetry provides limited protection 
from instabilities. We developed a random matrix model for the Hessian matrix 7{ at a generic 
critical point and obtained an analytic formula for its eigenvalue spectrum, finding that a signifi- 
cant fraction of the eigenvalues of "H are negative. Eigenvalue repulsion makes large fluctuations 
of the spectrum statistically costly, and by building on the theory of fluctuations of extreme 
eigenvalues — and through extensive simulations of the full Hessian matrix — we argued that 
the probability P of a large fluctuation rendering l-L positive definite is P oc exp(— ciV^), with 
p 1.5 and c ~ 0.3. 

Eigenvalue repulsion also controls the stability properties of approximately-supersymmetric 
critical points, at which the F-term F is small compared to the supersymmetric mass scale 
rrisusy In this regime, only the two eigenvalues corresponding to the Goldstino direction risk 
becoming tachyonic. We computed the two smallest eigenvalues to quadratic order in F/rrisnsy, 
and showed that mixing with the supersymmetric masses shifts these lowest eigenvalues to 
negative values. We then studied the probability of a fluctuation to positivity, through analysis 
of the corresponding univariate statistical distribution and through simulations of the full mass 
matrix. In the approximately-supersymmetric regime we found P oc exp(— cA^^), with p ^ 1.3 
and c ~ 0.1. 

We emphasize that the assumption that W and K are random functions — and in particular 
that their derivatives are independent random variables drawn from some statistical distribu- 
tion — is essential. There are, however, physically-motivated situations in which W and K are 
not general random functions of all of their arguments. An important example consists of two 
decoupled sectors: if Nh heavy scalars receive large supersymmetric masses, and supersymme- 
try is dynamically broken in a decoupled system of Nl lighter scalars at a much lower scale, 
then for a single vacuum configuration of the light fields the corresponding number of vacua 
of the full system can be exponential in N^- It seems reasonable to expect decoupling of this 
sort, into a 'degeneracy sector' at high scales, and a dynamical supersymmetry breaking sector 
at low scales, in a variety of compactifications. Restricting to the Nl light fields, our analysis 
suggests that the fraction of critical points that are metastable is proportional to exp(— cA"]^). 
For A^^ ^ Nl, the result is a mild reduction in the number of metastable vacua. 

Let us reiterate: our finding that an exponentially small fraction of critical points in a generic 
supergravity theory are metastable vacua in no way excludes the existence of a tremendously 
large landscape of vacua. There are two primary reasons, one conceptual and one quantitative. 



The conceptual reason is the possibility explained in §6.2| and reviewed above of a decoupled 
system (violating our assumptions onW, K) in which the vacuum degeneracy is ensured by Nh 
fields that receive large supersymmetric masses. The quantitative reason is that the values of 
c, p that we have obtained are not so large as to entirely overwhelm the vast number of critical 
points in flux compactifications. 

The methods and results of this work could be of use in understanding the statistical prop- 
erties of the moduli mass spectrum in string compactifications, and in guiding the search for 
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de Sitter vacua. One clear implication of our findings is that a direct search for explicit de 
Sitter vacua in systems with (9(10) or more fields and reasonably general W and K is likely 
to be frustrated by the appearance of tachyons. Correspondingly, the most promising regimes 
are those in which our assumptions are strongly violated, e.g. approximately-supersymmetric 
critical points for which the superpartner of the Goldstino enjoys special couplings to the re- 
maining fields. Understanding the incidence of such couplings in well-motivated supergravity 
theories, particularly those derived from string compactifications, is an important problem for 
the future. 
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In this appendix we briefly review some pertinent results from the the study of the distribution 
of non-supersymmetric vacua by Denef and Douglas P]. In particular, we will review how, for 
any fixed cosmological constant, the density of critical points — stable and unstable — grows 
linearly with F towards the boundary of the approximately-supersymmetric regime. This 
gives evidence for the expectation that a "generic critical point" typically does not exhibit any 
particular hierarchy between mgusy and F, i.e. typical critical points in random supergravity are 
not predominantly of the approximately supersymmetric kind. We also comment on how the 
distribution of metastable vacua is modified by the exponential suppression of the probability 
density found in ^ 

The density of critical points with cosmological constant {V) = v can be evaluated from 



Just as in reference j^j, we assume fiat prior probabilities for W, F, Z and U between and 
a cutoff A. Although this assumption is made here for simplicity, interesting domains of the 
string theory landscape have been argued to be well-described by these priors. The measure 
used in this appendix is 



A The distribution of critical points 




(A.l) 



dfx[W,F,Z, 



U]=C d^W d^^F d' 



(A.2) 
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where kz = N{N + 1), ku = fiN + 1){N + 2), and C is a normalization constant. 

As reviewed in ^ the critical point equation can be written as an eigenvalue equation 
for the matrix Ai of equation (2.5), enforcing that M has an eigenvalue equal to 2\W\ with 
F being proportional to the corresponding eigenvector. Thus, the critical point equation can 
be simplified by expressing the integration over Zab (and thereby Ai) as an integral over the 
ordered eigenvalues, Ai < . . . < Aat, and unitary rotations U. Let us denote an orthonormal 
eigenbasis of as e^, with corresponding eigenvalues ±Aa. In this basis F has components 

\F\r^t, with Ell {{vtYvt + iVaYVa 



1. In this notation 

N 



6''mM -2\W\)F 
1 



L n S{vn>^a-2\W\)) s{r^-i-\a-2\W\) 

I a=l 



N 



\2N 



a=l 



n 5k+(A.-2iiyi) 



Xa + 2\W\ 



The cosmological constant constraint 6{V — v) can be written as 



5{V -v) = -3|l^p 



Si\W\ 



w] 



3{\W\ 



w] 



(A.3) 



(A.4) 



where = | (F^ — v). Although the integral (A.l ) can also be estimated in the generic regime 
, this evaluation is slightly technical, and for the purpose of this appendix it 



in which F 



~ m 



susy 



suffices to discuss the approximately-supersymmetric regime of Denef and Douglas [2] . In this 

/N , and the determinant of the Hessian 



m. 



case, the integration over F is cut off before F 
appearing in the integrand of equation (A.l) is well-approximated by 

N 

n i^D" . 



susy 



IdefH 



2 2 



(A.5) 



a'=2 



with m\ as in equation (5.17). The number of critical points is 



N. 



crit.ptsA'^J 



n 



= C J dfil^F^F^UM] 

dXa6(v^{Xa-2\W\) 



a=l ^a-l 



d\W\^ 

\a + 2\W\ 



dF 5{\W\-w) 



F ?>{\W\+w) 
/(Ai,...,AAr) IdefHl , 



(A.6) 



where /(Ai, . . . , Aat) denotes the joint probabihty density of equation (3.8 ), e is a small number, 
and for notational convenience we have defined Aq = and Atv+i = A. In the approximately- 
supersymmetric regime, only Ai has a non- negligible probability density at 2|iy|. With this 
observation, the integral simplifies to 



A^crit.pts.(t;)=C /"rf/i[^^,(]^,[/,W](n<5^(r/+))(n5^(r^„-)) f 

a'=2 a=l "'0 



W 



X 



\i + 2w 



N 

n 



V+i 



d\n 



a' =2 V-1 



/(Ai = 2w, 



F 6w 
Xn) |detH|, (A.7) 
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where now the constant e' < 1 encodes the assumed hierarchy between Aa' and F . Since the 
probabihty density of the smallest eigenvalue exhibits a linear cleft for small arguments, cf. 



equation (3.9), we can (heuristically) write /(Ai = 2w, . . . , Aat) = 2kw /(A2, . . . , Atv), where k 
is an 0{1) constant. This simplifies the integral to 



crit.pts. \ ^ ) 



k c 



N 



N 



X 



e'\2 



dF 



2 2 
m , m_ 



n / 

,a'=2 -^V-i 



a'=2 



a=l 



A^) |det7/f/2^ 



(A.8) 



where H' denotes the truncation of "H to exclude the Goldstino direction. The scaling of the 
number of critical points with F at fixed cosmological constant is evidently determined by the 
factor 

[ ^ mlml. (A.9) 
Jo F 

For typical values of Um ~ rrisusy, the Goldstino masses-squared are each of order F, and 
the number of critical points scales with F as 



A^crit.pts.(^) ~ 



dF F. 



(A.IO) 



and thus, for any given scale of the supersymmetric masses, the critical points are more nu- 
merous towards the upper edge of the domain of approximate supersymmetry. 

We conclude with some simple remarks. This scaling of the number of critical points with 
F is consistent with the computation by Denef and Douglas of the scaling of metastable vacua 
with F, 



F^ dF, 



(A.ll) 



which is not surprising since the above computation closely mimics that of The different 
scalings of the number of critical points and the number of vacua can be understood from 
the additional fine-tuning necessary to obtain stability. In the approximately-supersymmetric 
regime it is necessary to tune |f/iii| ^ C>{F), which gives an additional factor of F"^ from 
the measure c?|f/iii| Furthermore, as reviewed in ^ the intent of this fine-tuning is to 

lower the scale of the Goldstino mass-squared to (9(-F^) in order to improve the probability of 



positivity of "H. By equation (A.9), this provides two more powers of F, from which equation 



(A.ll) follows. 



Finally, with these flat priors the additional A^-dependent (but F-independent) fine-tuning 
explored in this paper modifies the density of non-supersymmetric vacua in the approximately 
supersymmetric regime by 



lla'=2 Jo 



^ d\Una' 



la' 



Y[a'=2 lo 



N 



~ e 



d\Uua'\ \Ulla 



(A.12) 
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where iV", with a < 1, parameterizes the number of terms in S of equation (5.20) that need to 
be fine-tuned in order for a fluctuation to positivity of the smallest eigenvalue to become likely. 
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